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^-H Abstract 

o • 

ps| "A Generalized Master Equation (GME) is used to study the thermoelectric currents through a quantum dot in both the 
^ transient and steady-state regime. The two semi-infinite leads are kept at the same chemical potential but at different 
, ^ 'temperatures to produce a thermoelectric current which has a varying sign depending on the chemical potential. The 
Coulomb interaction between the electrons in the sample is included via the exact diagonalization method. We observe 
[a saw-teeth like profile of the current alternating with plateaus of almost zero current. Our calculations go beyond the 
I hnear response with respect to the temperature gradient, but are compatible with known results for the thermopower in 
the linear response regime. 
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1. Introduction 

] The electrical conduction of open nanoelectronic devices 
driven by electric potentials or fields generated in various 
ways is a major topic in mesoscopic physics. Outside this 
tarea complementary research on thermoelectric currents, 
thermopower, and related thermal properties in the quan- 
tum regime for systems like quantum dots has been more 
'difficult, but on a growing trend in the last two decades. 
Temperature control down to the milli Kelvin range and 
temperature gradients at nanoscale are attainable in the 
■laboratories and generate new scientific opportunities [1]. 
\ The thermopower of quantum dots was initially stud- 
ied theoretically by Beenakker and Staring [2]. They cal- 
culated the Seebeck coefficient 



S = — lim — — - 
AT^O AT 



(1) 



where V is the voltage generated across the quantum dot 
Weakly connected to electron reservoirs at a temperature 
bias AT, under the condition that the current between 
the two reservoirs is zero. They obtained oscillations of S 
around S* = as function of the Fermi energy in the reser- 
voirs, with symmetric positive and negative values. The 
Coulomb electron-electron interaction in the quantum dot 
was included in the charging ( "orthodox" ) model, and the 
result of it was a saw-tooth profile of the thermopower, 
with oscillations having the positive slope smaller than the 
negative slope. The predicted results were confirmed in 
subsequent experimental work by the same team [3] and 
also by Dzurak et al. [4]. Few years later Dzurak et al. 
published a new series of measurements which show that 
at temperatures below 100 mK the saw-tooth oscillations 
of the thermopower vs. the Fermi energy alternate with 
plateaus of zero thermopower [5]. A qualitative graph of 
the thermopower is shown in Fig. 1. 



While the saw-teeth were attributed to sequential tun- 
neling and high temperatures [2], the zero plateaus were 
initially attributed to many-body effects. Later on Turek 
and Matveev derived a theory of the thermopower of quan- 
tum dots in which the zero plateau at low temperature 
is a result of cotunneling [6]. A more complex cotunnel- 
ing theory, beyond the limit of weak tunnel coupling, and 
including quantum fluctuations, was proposed by Kubala 
and Konig [7], and later by Billings et al. who also included 
exchange effects [8]. Further experimental results were in- 
terpreted in terms of sequential-tunneling dominated ther- 
mopower at high temperatures, leading to a saw-tooth pro- 
file, and cotunneling onset at low temperatures, leading to 
zero plateaus [9] . In a recent experimental paper by Svens- 
son et al. [10] the thermopower of quantum dots is system- 
atically investigated and the lineshape is carefully analyzed 
in various conditions. Periodic sequences of a negative peak 
followed by positive peak followed by a zero plateau of the 
thermopower of the quantum dot as function of chemical 
potential in the leads are clearly seen over large intervals 
of chemical potentials. The interpretation of these results 
is done using a Landauer formula with an empirical trans- 
mission function. 



In the constant interaction model the effects of the 
Coulomb interaction are taken into account using only a 
charging energy, with no other effects on the energy spec- 
trum [11, 12]. The transport is described by a series of reso- 
nances, and in between the resonance peaks there are zero 
current plateaus. This is the standard Coulomb blockade 
picture [12]. In the Landauer-Biittiker approach this phe- 
nomena is described using a transmission function T{E) 
representing the resonances in the transport, usually a 
Lorentzian of width F centered on some resonance ener- 
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Figure 1: A qualitative behavior of the thermopower S and thermo- 
electric current / generated as a linear response to a temperature 
bias between leads, for a quantum dot, as function of the chemical 
potential in the leads for two regimes: high temperature, such that 
the resonance width F << ksT (upper panel) and low temperatures, 
i. e. F >> kgT . The red lines show the current and the green line 
shows the thermopower. Saw teeth alternating with zero plateaus are 
obtained for the current. The thermopower follows this behavior only 
at low temperatures. The chemical potential is in units of fe^T, but 
arbitrary units (a.u.) are used on the vertical axes. 

gies En, and the current is then given by [13] 

Hp) = I dET(E}(J(E,n} - f(E,n)) (2) 

where Tl ^ T ~ ^ and Tr = T + ^ are the tempera- 
tures of the two contacts (left and right) having identical 
chemical potentials fj,. There are two regimes that we want 
to consider from Eq. (3): (i) T < fcT and (ii) T > kT. 
In case (i) the current is non-zero around resonance peaks 
in an interval ~ kT. Between those intervals the current 
is exponentially suppressed due to the derivative of the 
Fermi function. In the other case, (ii), the current is pro- 
portional to the derivative of the Lorentzian peak, resulting 
in an suppression of the current between peaks. We 

assume both T and ksT smaller than the energy separation 
between resonances. 

The temperature driven current given in Eq. (3) will give 
rise to a voltage drop over the sample. In linear response 
the thermovoltage is given by 

VrM ^ ^, (4) 



where G is the conductance of the system. The conduc- 
tance in case (i) is determined by peaks of width ^ kT, 
with exponential suppression between adjacent resonances 
(assuming kT ^ AE) . The ratio of / and G will thus lead 
to a saw-tooth like pattern without plateaus in Vrh, even 
though there are plateaus in /. In case (ii) the conductance 
is proportional to the Lorentzian peak, which tends to zero 
like 1/E'^ outside the resonances. In this case the thermal 
voltage is saw-tooth like but with plateaus between reso- 
nances. 

In the present paper we use the generalized master equa- 
tion (GME) as a tool to understand the electric currents 
generated in an open quantum dot due to a temperature 
bias. The dot is connected to external leads seen as elec- 
tron reservoirs and kept at the same chemical potential. We 
obtain the currents in the leads produced by a finite tem- 
perature bias. We calculate the time dependent currents 
when the leads are gradually coupled to the sample and 
we find numerically the asymptotic currents in the leads 
in the steady state. To our knowledge the master equa- 
tion has not been commonly used for the thermoelectric 
response of open systems. One approach was performed by 
Koch et al. who used in fact rate equations, neglecting the 
off-diagonal elements of the reduced density operator, but 
including cotunneling [14]. Our method is not restricted to 
the linear response to a small temperature gradient. We ob- 
tain a line shape of the currents vs. the chemical potential 
as illustrated in Fig. 1 for low temperatures, where zero 
plateaus are also expected for the thermopower [5, 9, 10], 
as illustrated in Fig. 1. Our method however does not in- 
clude cotunneling effects, but only sequential tunneling. In- 
stead, the Coulomb interaction in the dot is completely in- 
corporated using the method of exact diagonalization. We 
also discuss time dependent and transient currents in the 
system and the effect of a third terminal attached to the 
quantum dot. The third terminal was proposed in order to 
create a phase-breaking mechanism inside the dot [15]. 

The paper is organized as follows: The model and the 
methodology are described in Section 2, analytical calcu- 
lations in Section 3, the numerical results are presented in 
Section 4, and the conclusions in Section 5. 



2. The Model 

The physical system consists of a sample connected to 
two leads acting as particle reservoirs. We shall adopt a 
tight-binding description of the system: the sample is a 
short quantum wire and the leads are ID and semi-infinite. 
The sample can also be seen as an elongated quantum 
dot. In this work we consider a sample of 3 sites. This 
number optimizes the computational time and the physical 
phenomenology which we intend to describe. A sketch is 
given in Fig. 2. The left lead (or the source, marked as L) 
is contacted at one end of the sample and the right lead 
(or the drain, marked as R) is contacted at the other end. 
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Figure 2: A sketch of the system under study. A ID lattice with 
3 sites ("the sample") is connected to two semi-infinite leads via 
tunneling. The left lead is connected to the left end of lattice, while 
right lead is connected to the right end. The contacts (xli Xr) ^^'^ 
time-dependent. 



The Hamiltonian of the coupled system reads as 

Hit) =. J2 He=L,R + Hs + Hrit) = Ho + Hrit) , (5) 

e 

where Hs is the Hamiltonian of the isolated sample, in- 
cluding the electron-electron interaction. 
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The (non-interacting) single-particle basis states have 
wave functions and discrete energies En- Hi, with 

{£} — {L, R), is the Hamiltonian corresponding to the left 
and the right leads. The last term in Eq. (5), Ht describes 
the time-dependent coupling between the single-particle 
basis states of the isolated sample and the states {ipqe} of 
the leads: 



(7) 



The function xe{t) describes the time-dependent switch- 
ing of the sample-lead contacts, while dj^ and Cqi cre- 
ate/annihilate electrons in the corresponding single- 
particle states of the sample or leads, respectively. The 
coupling coefficient 



VergimniH): 



(8) 



involves the two eigenfunctions evaluated at the contact 
sites (0, if), being the site of the lead £, and ig the 
site in the sample. The wave functions in the leads are 
'4'qi (0) = \/sm q/2T with r the hopping energy in the leads 
(the same for all leads), the energy spectrum of the leads 
being ef(q) — 2TC0sg [16]. The hopping energy in the sam- 
ple, denoted as ts, will be considered different than in the 
leads, and will be used as the energy unit. In Fig. 2 the 
left lead is connected to the site il — ^ and the right lead 
on the site i^j = 3. The parameter plays the role of a 
coupling constant between the sample and the leads. 

We will ignore the Coulomb effects in the leads, where we 
assume a high concentration of electrons and thus strong 
screening and fast particle rearrangements. The GME is 
formulated in the Fock space and therefore it is natural to 
include the Coulomb electron-electron interaction in the 
dot in a complete many-body manner and to calculate the 
exact many-body states using a straight forward diagonal- 
ization in the basis of occupation numbers. This method is 



known as exact diagonalization in the community of meso- 
scopic physics, but often called configuration interaction 
by quantum chemists. The many-electron states (MES) are 
calculated in the Fock space built on non-interacting single- 
particle states [17]. Since the sample is open the number 
of electrons is not fixed, but the Coulomb interaction con- 
serves the number of electrons, which means the interacting 
eigenvectors are linear combinations of the non-interacting 
eigenvectors with a fixed number of particles. 

The equation of motion for our system is the quantum 
Liouville equation. 



ihW{t) = [H{t),W(t)] 



(9) 



where W{t) is the statistical operator of the total system 
made by the sample plus leads, which are connected at 
time t = 0. Before the connection, at t < 0, the sample 
and the leads are independent and in equilibrium, meaning 
that W{t < 0) = plPrPs, i- e. simply the product of the 
density operator of the sample, ps, the left lead pL and the 
right lead pL- 

Following the Nakajima-Zwanzig technique [18] we de- 
fine the reduced density operator (RDO), p{t), by tracing 
out the degrees of freedom of the environment, the leads in 
our case, over the statistical operator of the entire system. 
Wit) 

p{t) = TrLTTBW{t). (10) 

For a sufficiently weak coupling strength (Vt) one obtains 
the non-Markovian integro-differential master equation for 
the RDO 



m 



Pit)] 



where the operators Qqi and Hqg are defined as 



(11a) 



Uq.is) = e-^- (r>(s)(l - /,) - pis)T^Ji 



(lib) 



and /£ is the Fermi function of the lead i describing the 
state of the lead before being coupled to the sample. The 
operators Tqe and T^^ describe the 'transitions' between 
two many-electron states (MES) when one electron enters 
the sample or leaves it. 

The GME is solved numerically by calculating the ma- 
trix elements of the RDO in the basis of the interacting 
MES, in small time steps, following a Crank-Nicolson al- 
gorithm. More details of the derivation of the GME can 
be found in Ref. [16]. The calculation of the interacting 
MES is described in Ref. [17]. The switching functions xt{t) 
must be defined. For example any function starting at zero 
and gradually increasing to one can be used to obtain the 
steady state in the asymptotic limit. In principle any other 
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time dependence can be used, like steps or periodic func- 
tions. 

Mean values of observables can by obtained by taking 
the trace of product of the corresponding operator and the 
RDO. The total time dependent charge in the sample is 
found by using the number operator M = <ij„<im: 

{Qit))=eTT{pAf} = eJ2Nj2{aN\p{t)\aN), (12) 

N an 

where ajv denotes the (Coulomb interacting) MESs with 
fixed number of electrons N . Remark that one can also cal- 
culate the partial charge accumulated on A^-particle MESs. 

The currents in the leads are then found by taking the 
derivative of Eq. (12) with respect to time. The time deriva- 
tive of the RDO can be substituted by the right-hand side 
of the GME [Eq. (11a)] and so it is possible identify the 
currents in each lead, 



N QJV 



+ h.c. 



(13) 



3. GME With One Energy Level 

To better understand the thermal effects in the GME 
we solve the GME for a sample with one single site (a 
single level quantum dot). To simplify Eq. (11a) we start 
by inserting 



p{s)^Ul{t~s)p{t)Usit-s) 



(14) 



i.e. we propagate the density matrix backwards in time 
and take it outside the time integral in Eq. (lib). In fact 
this is the Markov approximation. 

The Fock space contains now only two states, the vac- 
uum state with energy and a single particle states with 
energy E > 0. We use here E = 2 units of tg. Of course, for 
a one-site model the hopping energy has no meaning, the 
result being only the diagonal term in the lattice Hamil- 
tonian. We want to find the occupation of these states, 
i. e. the diagonal elements of the RDO, for this two-level 
model. In this case Eqs.(ll) are greatly simplified because 
only one of the transfer matrix elements (Tqi) is non-zero: 



(i|r,,|o) = ^rJ„(i|4lo)=T,V 



(15) 



This simplifies the commutator Eq. (11a) and gives 

:{0m^-^J2 J dqRc{m,,\l){l\T,i\0)}. (16) 



The trace of the RDO is 1, therefore the other diagonal 
element is given by (l|/o|l) = 1 — (0|p|0). Only one matrix 



element for fl is needed to find the diagonal elements of 
the RDO 



(o|rJ,|i) 



: {e-" 



-it{et{q)-E) 



1 



eiiq) - E 

X ((1H1)( 1 - /,(e,(<z)) ) - (0|p|0)/,(e,(g)) 



(17) 



The g-integral can be evaluated in the steady-state limit if 
we admit that for t oo the real part of (0|f2g£|l) contains 
a delta function S {^i{q) — E^ . With this approximation we 
actually neglect the level broadening F due to the lead-dot 
coupling [19]. That broadening will be discussed in the 
next section. In the steady-state limit the diagonal values 
of p approach a constant value and thus the right hand 
side of Eq. (16) must be zero. Using Eq. (8) and the delta 
function one obtains 

^ [(l|p|l) (1 - ME)) ^ {0\p\)ME)] = . (18) 



Therefore the occupations in the steady-state are 



(0|p|0) -1- 
and the current is 



(l|p|l) = l-(0|p|0), (19) 



Ji 



^{{0\p\0)ME) - (l|p|l)(l - ME))}. (20) 



The equation for the occupation and the currents in the 
leads was derived using the fact that in the steady-state p 
is constant and thus its time derivative is zero. This implies 
that in the steady state the currents in the left and right 
leads are equal: what goes into the sample must also go out. 
Therefore using Eq. (19) we can eliminate p from Eq. (20) 
and we obtain: 



Jl 



1 



1/2 



(ME) - fn{E) 



fR{E)-fL{E) 



(21) 



Which makes it clear that in the steady-state the currents 
are equal, but with opposite signs Jl = ^Jr, with positive 
"in" from the left led, and negative "out" into the right lead. 

The results of the analytical calculations can be seen 
in Fig. 3. In Fig. 3(a) we show the results for a thermal 
bias, i. e. for different temperatures in the leads, Tl > Tfj 
and PL = PR, and also for an electrochemical bias, i. e. 
Tl = Tr and pL > PR- Fig. 3(b) shows the charge in the 
system calculated using Eq. (19). The current generated 
thermically can be positive or negative, the transition oc- 
curring when the Fermi level is equal to the energy of the 
single level Pj^ = E. At this point the state is half filled 
and the current in the leads is zero. The probability of a 
transition to/from the energy level to/from the left and 
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Figure 3: Analytical calculations for a system with one site, (a) The 
red solid line shows the currents in the left lead for a thermal bias, but 
no electrochemical bias. The blue dashed line shows the current for an 
electrochemical bias, but no thermal bias. The thermal bias used in 
the calculations was ksT^ = 0.25, kgTn = 0.10, with fj, = fij^ = fiR. 
The electrochemical bias was A/i = 0.15, with kgT]^ = kgTji = 0.10. 
(b) Charge in the system in the two cases. On the horizontal axis 
we use the chemical potential fi = fiL = MH fo'' the thermal bias and 
the mean value fi = (/i^ + iiji)/2 for the electrochemical bias. The 
green vertical line shows the energy of the resonance, E = 2. 

right lead is equal. When the state is less than half filled 
a transition from the energy level to the right lead has a 
higher probability than to the left lead. This gives a pos- 
itive current flowing from the left to the right. Once the 
state is more than half filled a transition from it to the left 
lead is more probable. Giving a negative current flow from 
the right to the left. This difference in transition probabil- 
ity is due to different temperatures at the left and right 
lead. 

According to Eq. (21) the current is essentially the differ- 
ence of two Fermi functions. For the thermoelectric current 
they are centered at the same chemical potential, but they 
have different widths. Tor the pure electric current they are 
centered at two different chemical potentials, but they have 
the same width. Therefore the width of the current peaks 
in Fig. 3 are only due to the temperatures or chemical 
potentials. For example the width of the current created 
by the chemical potential bias, at half height, should be 
A/.t -I- 2kBT = 0.15 -I- 2 X 0.1 = 0.35, which is consistent 
with Fig. 3. (The level width F has been neglected.) 

We can conclude this section with the idea that the ther- 
moelectric currents can be understood as being related to 
the difference between two Fermi functions with the same 
center but different widths. Therefore, as pointed out re- 
cently by Tagani and Soleimani [20] , there are two reasons 



for the current to become zero: (1) half filling, where the 
two Fermi functions are equal to 0.5, and (2) integer fill- 
ing, where both are or 1. This idea will help us to under- 
stand more complex results incorporating many-body ef- 
fects. For comparison the currents generated by an electro- 
chemical potential bias are given by the difference between 
two Fermi functions of the same width but with different 
centers, and thus have a constant sign and a maximum at 
half filling. 

4. Results for a Many-Body System 
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Figure 4: Electrochemical potential diagram for a system with 3 sites 
and Coulomb interaction strength Uc = 2.0. Single particle states 
blue dots ( — • — ). Two particle states red squares ( — ■ — ). Three par- 
ticle state brown diamond (—♦—). 

The MESs of the sample are characterized by the elec- 
trochemical potentials /i^-* :— s'"^ — f^^j^, where f^'' is 
an energy of the sample spectrum containing N parti- 
cles, i — indicating the ground state and i > the ex- 
cited states. In Fig. 4 we see the electrochemical po- 
tential diagram for a system with three lattice points 
and Coulomb interaction strength Uc = 2.0 (units of tg). 
The system has three single particle states with energies 
El = 2- y(2) w 0.59, E\ = 2,Ef = 2 + V2~ 3.41, three 
two-particle states with energies £'2 = 4, i?2 = 6, i?| = 7 
and one three-particle state with energy E^ — 11. For 
the single-particle states (iV = 1) the chemical potentials 
are in fact the single-particle energies. The effects of the 
Coulomb interaction on the two (N — 2) and three-particle 
states {N = 3) is to shift them upwards. 

Fig. 5 shows the current in the left lead in the steady 
state for both Uc — 0.0 and Uc — 2.0. The green vertical 
lines represent the chemical potentials in Fig. 4. The blue 
line (--•--) shows the current for a system with no Coulomb 
interaction. The current is zero at the points fi « 0.5, ii — 2 
and n « 3.5. These points correspond to half filling of 
states. In between such points the current is again zero 
when integer filling occurs (e. g. for fi w 1.2). The fill- 
ing can be seen in Fig. 6(a), were the charging of the sin- 
gle, two- and three-particle states is shown for the non- 
interacting case. At the point « 0.5 the single particle 
states are half charged, and at /i = 2 the two-particle states 
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Figure 5: Comparison of the current in the left lead with and with- 
out Coulomb interaction. Red solid line ( — ■ — ) has the Coulomb 
strength Uc = 2.0 while the blue dashed line (-■•--) has Uc = 0.0. 
The temperature is kgTi^ = 0.25 in the left lead and kgTn = 0.10 
in the right lead. The lead-dot coupling parameters are Vi = Vl = 
Vn = 0.75 and the hopping energy in the leads is t = 4 (units of ts). 



are half charged and the single particle half discharged. 
The last point /i « 3.5 is the point is where the two-particle 
states are half discharged and the three-particle state half 
filled. The red line (— ^) shows the current for the same 
system with Coulomb interaction {uc — 2.0). The effects of 
the Coulomb is to create plateaus of zero current at /x = 2 
and /i « 5.5. These points correspond to integer filling of 
states. 




Figure 6: Charge as a function of the chemical potential. Brown solid 
line ( — ■— ) shows charging for the single particle states, black dot- 
ted lines (■■■A-- ) for two-particle states, violet dashed lines (-■•-) 
for the three-particle state, and blue dashed lines (--•--) show 
the total charge, (a) Without Coulomb interaction, fcsT^ = 0.25 
and fcsTi = 0.10. (b) With Coulomb interaction {uc = 2.0), 
kBT^ = 0.25 and fc^T/, = 0.10. (c) With Coulomb interaction (lie = 
2.0), kBTL = 0.55 and kgTL = 0.40. 



Fig. 6(b) shows the charging for the system with 
Coulomb interaction {uc — 2.0). The blue line (--^-) shows 
the total charge and has integer fillings at /j, = 2 and ^ = 5. 
At ^ = 2 the single particle states (— ■— ) are completely 
filled and at = 5 the two-particle states (■■■a- ) are filled. 
While the points /i « 0.5, /i sa 3.5 and n — 7 correspond 
to half filling. 

In Fig. 7 we show the current for two different tem- 
peratures. The red dashed line (-■•--) has the tempera- 
tures ksTL = 0.25 and fcsTi = 0.10 in the left and right 
leads respectively. The blue solid line (— ■— ) has the same 
temperature difference between the leads A/c^T 0.15 
but the temperature has been increased by 0.35 in both 
leads. The temperature in the left lead is then /cbTl — 0.55 
and ksTji — 0.40 in the right lead. One effect of increasing 
the temperature in both leads is that the current increases 
in magnitude. The zeroes in the current due to half filling 
are not affected but the plateaus due to integer filling are 
raised or lowered and reduced. This behavior can be also 
explained using the charging diagram. Comparing Fig. 6(c) 
with Fig. 6(b) one can see that around chemical potential 
/X = 2, due to the increased temperature, the population 
of single particle states decreases considerably from one, 
but the population of the two-particle states increases con- 
siderably from zero. Therefore the single particle states 
create a negative current whereas the two-particle states 
create a positive one (as illustrated in Fig. 3). The sum 
of these two contributions create the total current, in this 
case positive. For chemical potentials around 5.5 the total 
current is a combination of two-particle and three-particle 
components. The population of the two-particle states is 
less than one, but larger than one half, hence with a nega- 
tive contribution to the current, whereas the three-particle 
population is between zero and half, hence with a positive 
contribution. In this case the total current is negative. 
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Figure 7: A comparison of the current in the left lead for two different 
temperatures. The red dashed line (-■•--) has kBT^^n = (0.25, 0.10). 
While the blue solid line (-■-) has fcsTi^fl = (0.55, 0.40). The 
temperature difference between the leads is the same in both 
cases AfcsT = 0.15. 

It is interesting to compare the current created by a 
temperature bias with that due to an electrochemical bias 
in the many-body case. This is what we show in Fig. 8. 
The current driven by the thermal bias is the same as in 
Fig. 7. The temperature and chemical potential differences 
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Figure 8: A comparison of the current created by a temperature 
bias with the current created by a chemical potential bias. The 
red solid line (— • — ) shows the thermal current with ksT^ = 0.25, 
kgTn = 0.10 vs. the chemical potential fi = /x^ = /ij^. The blue 
dashed line (-■•-) shows the electrical current with no tempera- 
ture bias, kgTi^ = kgTn = 0.10, but Afi = fij^ — fiji = 0.15, i. e. en- 
ergetically identical to the the previous thermal bias. On the hor- 
izontal axis we use the mean value of the left and right chemical 
potentials /i = (/x/^ -|- /ifl)/2. The lead-dot coupling parameters are 
Vi = Vl = Vr = 0.75 and the hopping energy in the leads is r = 4 
(units of ts). 

between the left and right leads are identical on the en- 
ergy scale. Obviously in the later case the current is pos- 
itive, with peaks indicating the half filling of the ground 
state with N electrons, and with zero values in between 
these peaks, indicating the Coulomb blocking of the trans- 
port. We also obtain a small peak at /i = 2 which corre- 
sponds to a small current going through the first excited 
single-particle state. The magnitude of the Coulomb inter- 
action is incorporated in the energetic separation between 
the main peaks and may be estimated from the energetic 
length of the zero plateau of the thermoelectric current. 

Now we can estimate the level broadening due to the 
coupling between the leads and the dot, F. For exam- 
ple, the first peak of the current produced by the chem- 
ical potential bias around /i — 0.6 has a width at half 
height of 0.58 (units of ts). As shown in the previous sec- 
tion the broadening corresponding to the bias window is 
-I- 2kBT = 0.35. The difference can be attributed to 
the coupling, i. e. F = 0.58 - 0.35 = 0.23. This results 
is consistent with calculations at lower bias windows and 
temperatures (not shown) . In principle the broadening due 
to coupling is proportional to the coupling coefficients, Eq. 
8, which are not simple parameters, but a matrix elements 
and functions of energy which are hidden in the numerical 
solution of the GME. Nevertheless, in general, from previ- 
ous studies [16, 17, 21], we expect F a fraction of Vi, like 
25-50% or so. Comparing the estimated F — 0.23 with the 
thermal energies used 0.10 < fc^T < 0.55 we can conclude 
that our results correspond to an intermediate parameter 
regime, in between the situations presented in Fig. 1. 

Another comparison we want to make is between the 
complete Coulomb effects, which we describe via the exact 
diagonalization of the sample Hamiltonian, and the charg- 
ing model which is usually invoked in problems related to 
the Coulomb blockade. The latter is known as the "ortho- 



Figure 9: The blue solid line ( — ■ — ) shows the current in the left 
lead calculated using the GME with Uc = 2.0. While the red dashed 
line (-*--) shows the current calculated with the GME using eigen- 
vectors corresponding to non-interacting electrons, but with an en- 
ergy spectrum forced to be that of the system with Uc = 2.0. 

dox" model and it assumes that the energy of a two-particle 
state is the sum of the energy of the first two single particle 
states plus some charging energy estimated as the electron 
charge divided by the capacitance of the quantum dot [12]. 
This is essentially a mean field assumption. In order to 
compare the two approaches we calculate the current gener- 
ated by the temperature bias by using the non-interacting 
eigenvectors in the many-body sample Hamiltonian (i. e. 
the basis in the Fock space), and the ground state energy of 
the sample for 2 electrons as given by the exact treatment. 
With this ansatz we describe uncorrelated particles with 
the energies of the exact states. The results are shown in 
Fig. 9. In the "orthodox" model the plateau of the thermo- 
electric current, at least for a number of electrons between 
one and two, softens or vanishes, possibly depending on 
strength of the Coulomb interaction. Our interpretation of 
this result is that the zero plateaus of the thermoelectric 
currents are at least partially and effect of electron-electron 
correlations, not captured by the "orthodox" model of the 
Coulomb blockade. 

Next we show results for a sample attached to three 
leads, i. e. with a third lead connected to site number 2. 
Such a setup was proposed recently [15] in a different con- 
text, where the third terminal was used as a phase breaking 
mechanism of the ballistic electron propagation between 
the other two terminals. In our case the third lead has 
a different effect: it softens the Coulomb blocking of the 
currents by allowing the electrons entering the sample to 
diffuse into the third terminal. In Fig. 10 one can see how 
the zero plateaus of the thermoelectric current vanish in 
the presence of the third terminal, indicating the additional 
degree of freedom offered to the electrons injected into the 
sample. 

Finally we also show time dependent currents generated 
by a temperature bias. Until now we discussed only the 
currents in the steady state obtained by time integration, 
i. e. in the asymptotic time limit. But the important ad- 
vantage of the GME is that one can examine the transient 
regime in the system or another kind of time evolution. 
Fig. 11 shows time dependent currents of a turnstile model 
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Figure 10: A comparison of the current in the leads between systems 
with two and three leads. The blue dashed line (--•--) shows the 
current for a system with two leads, as before, with ksTh = 0.25 
and kgTfi = 0.10, and the same chemical potential /i. The red solid 
line ( — ■ — ) shows the current in the same system, but with a third 
lead attached to the middle point in the sample with KgT]^ = 0.10 
and the same chemical potential as the other leads. 
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Figure 11: The current in the left lead for a turnstile system. 
The switching function is x^(*) = sin(ajt + 4>i) seen at the bot- 
tom, (a) Electrochemical bias /x^ = 3.25. (b) Electrochem- 
ical bias liL,R = 3.75. Other parameters are Uc = 2.0 and 
kBT^L.R) = (6.25,0.10). 



with a thermal bias. The temperature of the left lead is 
keTL = 0.25 and of the right lead fesT^ = 0.10. The 
switching functions Xii^) start at zero and continues with 
sinusoidal pulses, as shown n the lower panels of Fig. 11. 
A phase shift of half a cycle between the left and right 
lead is included. The current is in both leads are shown 
for two cases, hl = jJ-R = 3.25 and = IJ-r = 3.75. Ear- 
lier, in Eq. (21), we used the sign rule such currents which 
enter into the sample or exit from it have opposite signs. 
In this example we prefer to consider positive the currents 
from left to right and negative those from right to left. 
The currents follow the shape of the pulses applied to the 
contacts. Initially the pulses are large because the system 
starts out empty. The initial charging of the system is from 
both leads, with a positive current from L into the sample 
and a negative current from R also into the sample. After 
the transient phase the pulses stabilizes and the system en- 
ters in a periodic (steady) state. In Fig. 11(a) the current 
pumped in each cycle is positive whereas in Fig. 11(b) the 
current pumped over each cycle is negative. This is quali- 
tatively different from an analog pumping in the presence 
of chemical potential bias with the same temperature in 
the leads. In that case only positive net current can be 
transported over the sample in one complete cycle [21]. 



erated by a temperature bias between two leads in qualita- 
tive agreement with the experimental results for the ther- 
mopower of quantum dots at low temperatures [10], with 
saw teeth alternating with zero plateaus. The plateaus of 
the thermopower have been explained by other authors 
with the cotunneling processes across the leads through 
the dot [6] . Our formulation of the GME does not include 
cotunneling, but only the sequential tunneling. Therefore 
in our model the zero plateaus of the current can be ex- 
plained using the level spacing, the Coulomb blocking, and 
electron-electron correlations. Cotunneling effects may in- 
deed contribute to the plateaus of the current and includ- 
ing them into the GME for a multilevel system is a com- 
putational challenge and the goal of a future work. 

Our study may also be compared with a recent work on 
the thermoelectric properties of serially coupled quantum 
dots [20]. Those authors used rate equations and a short- 
range Hubbard model for the electron-electron interaction 
and obtained the saw teeth of the thermopower. Our sam- 
ple is designed as a series of sites, and therefore each site 
can also be interpreted as a single quantum dot with a 
single bound state. Therefore our model can be naturally 
used for serially coupled dots. 



5. Concluding remarks 

In this work we performed an exploratory study of 
the current driven through an open quantum dot in a 
thermopower setup using the generalized master equation 
(GME). The quantum dot has been connected to two leads 
having different temperatures and the same chemical po- 
tential. The quantum dot has been defined with a finite 
number of discrete states. The GME method allows us to 
describe time dependent, transient, and also steady states, 
by numerical integration in time. We also included the 
electron-electron interaction in the sample using the ex- 
act diagonalization procedure. We obtained currents gen- 
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